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^ Abstract 

,_i The acceleration of polarized electron beams in the blowout regime of plasma-based acceleration 

Qh is explored. An analytical model for the spin precession of single beam electrons, and depolariza- 



Q^ 



tion rates of zero emittance electron beams, is derived. The role of finite emittance is examined 



O 
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j^ numerically by solving the equations for the spin precession with a spin tracking algorithm. The 

O 

■^ analytical model is in very good agreement with the results from 3D particle-in-cell simulations 

>-. 

<~| in the limits of validity of our theory. Our work shows that the beam depolarization is lower 



for high-energy accelerator stages, and that under the appropriate conditions, the depolarization 
associated with the acceleration of 100-500 GeV electrons can be kept below 0.1 — 0.2%. 
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I. INTRODUCTION 

Polarized particle beams are widely used in state-of-the-art high-energy physics (HEP) 
experiments [IH3] . The beam polarization, which corresponds to the average of the particle 
beam spins, is an essential beam feature, critical to several fundamental physical problems, 
including precision tests to the standard model |3], the search for the Higgs boson [1], 
and the discovery of new physics, for instance, in higher dimensions |3]. State-of-the-art 
and future accelerators are aiming at the acceleration of highly polarized electron beams. 
For instance, design for the International Linear Collider (ILC) considers 0.5 TeV electron 
beams with polarizations above 80 % (close to the maximum polarization of 100%, where 
all the beam particles spins are aligned in the same direction), with typical depolarizations 
below 0.1 — 0.2 % [5], which are mainly associated with the beam-beam interaction at the 
interaction point. 

The maximum acceleration gradients that can be attained in standard accelerators, and 
that directly impacts on their size, and cost, is limited by material breakdown thresholds. 
The use of plasma waves to accelerate particles can, therefore, play an important role in 
future generations of accelerators by sustaining acceleration gradients which can be more 
than three orders of magnitude higher than standard accelerators [6] . 

In plasma based acceleration (PBA), accelerating structures can be excited by either 
particle or laser drivers. The most striking experiments using electron beam drivers [plasma 
Wakefield acceleration (PWFA)] doubled the energy of 42 GeV electron beams in 80 cm long 
plasmas [7j. These results were obtained in the blowout regime, where plasma electrons are 
evacuated from the region where the driver propagates. The resulting wakefield structures 
are characterized by linear accelerating and focusing forces, and are ideally suited for electron 
acceleration [SHTO]. 

The most important advances using laser pulse drivers [laser wakefield acceleration 
(LWFA)], also occurred in the blowout regime, and demonstrated the acceleration of 1 GeV 
electron beams in cm scale plasmas pTI - ITS] . Several techniques were also devised in order 
to control the acceleration processes [IE], adjusting the bunch energy, charge [T7H20] . and 
transverse features [21]. 

Recently, the use of proton drivers was also proposed as a means to accelerate electron 
beams even further to > 0.5 TeV in < 1 km plasmas [22j. These results show, clearly. 



the capability of PBA, and in particular, the capability of the blowout regime of PBA, to 
accelerate high-quality electron beams to high energies p3| ^^ in controlled acceleration 
scenarios. However, for a future plasma-based linear collider, or for a plasma afterburner [3, 
[251428] . it is important not only to achieve high final output energies, but also to achieve 
low and controlled depolarization rates. 

The depolarization of an electron beam is fully determined by the temporal evolution 
of the direction of the spin, s, of each beam electron. During its acceleration, the spin of 
the electron can change according to two different mechanisms: the stochastic spin diffusion 
from photon emission [2^ I5U] , and the spin precession around the electric and the magnetic 
fields [21 [31]. On one hand, the stochastic spin diffusion is a non-deterministic process 
consisting in the rotation of the spin vector in the presence of a magnetic field, with the 
emission of a photon. The Sokolov-Ternov effect p9] (spin-flip) is a well-known stochastic 
spin diffusion mechanism, which can polarize electron beams in circular accelerators. On the 
other hand, the spin precession is a deterministic process, and can be examined by treating 
the spin as an intrinsic electron magnetic moment. 

The spin precession is an important mechanism that has a decisive role in the design 
of accelerators. In circular accelerators, for instance, the spin precession can depolarize 
the beam completely if the spin precession frequency is a multiple of the orbital frequency 
(depolarizing resonance) [H [2]. In linear accelerators (linacs) and colliders, however, the 
typical depolarizations are as low as 0.1 — 0.5% for 0.5 — 1 TeV accelerators and colliders [3]. 
This work is then focused on the dynamics of spin precession in plasma wakes, while the 
nondeterministic spin diffusion is left for a future work. 

This Paper shows that the electron beam depolarization during the acceleration in PBAs 
may fulfill the requirements for high energy physics experiments, provided that the beam 
emittances are sufficiently low, and shows that depolarizations of 0.1 — 0.2% can be achieved 
in 100-500 GeV accelerators. In Section [III the Thomas-Bargman-Michel-Telegdi (T-BMT) 
equations, which describe the spin precession of relativistic charged particles, are used to 
derive a set of coupled equations for the spin precession in plasma accelerators. This analysis 
reveals that the spin precession is fully determined by the transverse forces that act on 



relativistic particles. In Sec. Ill, analytical expressions for the depolarization associated 



with zero-emittance electron beams are derived. In Section IV it is shown that externally 



guided propagation regimes lead to lower depolarization rates in comparison to self-guided 



propagation regimes. In addition, it is found that the depolarization of higher energy plasma 
accelerators is lower than that of lower energy plasma accelerators. In Sec. |V] the model 
is compared with numerical simulations in scenarios that go beyond the validity limits of 
the analytical theory. Two approaches were used: in the first approach, the prescribed 
electromagnetic fields and trajectories of electrons in the blowout regime are considered; the 
second approach uses the fully self-consistent electromagnetic fields and electron trajectories 



from 3D PIC simulations in QuickPIC [32]. Finally, in Sec. VI the conclusions are stated. 



II. SINGLE ELECTRON SPIN DYNAMICS 

In order to determine the polarization of the beam, the spin precession of a single electron 
is firstly examined. For non-relativistic electrons, the spin precession follows: 

ds s X B . . 

where /x = —g s/2 is the electron intrinsic magnetic moment, g ~ 2.0023228 is the dimen- 
sionless magnetic moment of the electron (g-factor), t is the time normalized to the inverse 
of the plasma frequency cOp = A/47re%o/m^, with e and rUe being the electron charge and 
mass, uq the plasma density, B is the magnetic field normalized to nieCUp/e, and where c is 
the speed of the light. 

Eq. ([l]) indicates that the spin of a non-relativistic electron - with relativistic gamma 
factor 7 = l/-\/l — v^ > 1, with v denoting its velocity normalized to c - precesses around 
the magnetic field lines with frequency fl^^ = — /xB = gB/2. However, most interesting 
particle acceleration scenarios use ultra-relativistic electron beams with 7^1. In order to 
obtain a correct description of the spin precession dynamics for this case, the relativistic 
generalization of Eq. ([I]), which is not covariant, is required, and given by 



ds 

dt 



a + -)(B-vxE)- v^^v ■ B 
7/ 7 + 1 



X s = f2 X s. (2) 



All the quantities in Eq. (pi) are in the laboratory frame, except for s, which is an intrinsic 
property of the electron and is, for that reason, described in the electron rest frame. More- 
over, the precession frequency fi is $7 = [{a + 1/7)(B — v x E — va7v ■ B/(7 + 1)]. The 
quantity a = {g — 2)/2 ~ 0.0011614 is the anomalous magnetic moment of the electron, and 
E is the electric field normalized to meOJp/e. 



Although the T-BMT equations are strictly valid when E and B are homogenous, they 
can still be used as long as Stern-Gerlach type forces(F^~*^) can be neglected P], i.e. as long 
as the spin dynamics does not change the electron trajectories. These forces are proportional 
to F^^'^ oc V(/i ■ B) oc {gh/2)rQijJi3B, where B oc rocos(co'/3t) was assumed, and where w^? = 
bjpj ^/^ is the betatron frequency. Thus, Lorentz forces, proportional to F^ oc E + v x B oc 
ro ^ i^s-G ^ Yir^j ^^ dominate the electron dynamics in the plasma wave, validating 
Eq. g for PBAs. 

Investigation of Eq. ([2]) in cylindrical coordinates provides a clear understanding of the 
physics of the spin precession for a single electron. In cylindrical coordinates, the spin vector 
is s = (sr, s<^, S2), where Sj = s ■ ei, and where ei correspond to the unit vectors in the radial 
(e^), azymuthal (e,^), and longitudinal (e^) directions. In addition, (r, 0, z) are the radial, 
azimuthal and longitudinal coordinates. To simplify the notation, each spin component is 
normalized to the absolute value of the electron spin, such that s takes values between -1 
and 1. 

The electric field inside the plasma wave has both radial and longitudinal components 
E = ErG,r + EzBz, with Ez ~ Er- Morcovcr, for electrons with 7^1, Vr '^ Vz — 1, which 
means that the first term on the right-hand-side of Eq. ^ coincides with the radial plasma 
focusing force Fr felt by a relativistic electron: 

(B - V X E) = (S^ - VzEr + VrEz) 6^ ~ {B^ - Er) 6^ = FrS^, (3) 

where Vr and Vz are the radial and longitudinal components of the electron velocity, and —F^ 
is the radial force that acts on a relativistic electron in a plasma accelerator. Furthermore, 
the magnetic field of cylindrically symmetric plasma waves is purely azimuthal, i.e. B = 
B(f,e^ [8]. Since for ultra-relativistic electrons 7^1 (the most relevant scenario for PBAs), 
V(f) "t^ Vz, where f^ is the azimuthal component of the electron velocity, then v ■ B = v^Bff, ^ 
1. Consequently, the last term on the right-hand-side of Eq. ^ can be neglected, and hence 
Fr fully determines the spin precession of a single electron in PBAs. 

Inserting Eq. dsl) in Eq. ([2]), and neglecting the v ■ B term in Eq. (|2| yields: 



dSr f 1 ^ 

= s^ a + 
dt V 7y 


Fr + S 


ds, / n 

dt V 7^ 


FrSr 


dt 





(4a) 
(4b) 
(4c) 

where = d0/dt. 

Eqs. Q represent a set of coupled differential equations which describe the spin precession 
of a single electron in the field structure of a PBA, valid as long as plasma waves are 
cylindrically symmetric. According to Eq. Q, the spin precession is mostly due to the 
anomalous electron magnetic moment when 7 ^ 1/a ~ 10^. For 7 ^ 1/a, the anomalous 
electron magnetic moment can be neglected. 

Solutions for Eqs. Q can be retrieved by calculating the transverse electron trajectory 
in linear focusing and accelerating forces. Linear focusing forces are ideal for particle ac- 
celeration, as they preserve the transverse emittance of the accelerated beam. This is very 
important for HEP experiments which require low emittance beams in order to achieve high 
luminosities. Both the blowout, weakly- relativistic and even linear-regimes are character- 
ized by linear focusing forces as long as the electron beam propagates close to the axis. 
Thus, Fr = F_L ■ ir = axx ■ xi^/|x_l|, where F_l = {Fx,Fy) and x_l = {x,y). The parameter 
a = 1/2 if the blowout is complete. Otherwise, a < 1/2. 

In the presence of uniform accelerating fields, -Eacceh the beam relativistic factor varies 
with 7 = 7o + -Eaccci^, where 70 = •yit = 0). This simple model for the acceleration of 
electrons in plasma waves is accurate for a significant fraction of the acceleration distance, 
and only fails when -Eaccei is significantly reduced due to the laser pulse pump depletion in 
LWFAs, or electron beam head-erosion in PWFAs, and when the acceleration distance is 
comparable with the dephasing length in both LWFAs and PWFAs. 

Under the conditions mentioned above, the radial electron oscillations in the plasma wave 
ion channel are given by 



(70 + ^acceit)i± (t) + 7oi± (t) + axx (t) = 0, (5) 

where x_lo = {xo,yo) is the initial transverse position of the electron. Eq. ([s]) can be solved 
analytically for an electron initially located at x^Q) with initial transverse momentum p^g = 



7v_Lo, for which ^/orfy/E^ccc\ > 1 and ^/ori / E.^ccc\ > 1, yielding: 



\l{t) J y -Spaced -Spaced J 

(a^77o) \ -2'accel -traced / 

Provided that x^o and p_Lo are known, Eqs. Q and ^ fully describe the spin precession 
dynamics of a single electron in a PBA. 

The most interesting limit for HEP applications is associated with low beam emittances, 
characterized by low average |p±o| iu comparison to the average |x_|_o|, such that the eccen- 
tricity of the electrons trajectory in the plane perpendicular to the propagation direction is 
much higher than one. According to Eq. ([6| this occurs for (|x_lo|) ^ {\p±o\)/ y/^'jo, where 
(|x±o|) — o"r is the typical transverse width of the electron beam, and where (|p±o|) is the 
typical electron beam perpendicular momentum. Equivalently, the effects due to finite beam 
emittance can be neglected if the normalized beam emittance, e^, is much lower than 

ernmm ■ mrad] < 188^^ ( ^^] ' ( ^' ] ( .!''' ^ '^' , (7) 

where the normalized emittance is estimated according to eAr = (|x^|)(|p_l|), and where Eq is 
the initial beam energy [we note that eN = Cn*^^ corresponds to matched PWFA propagation 
regimes which minimize the betatron oscillations of the driving electron bunch]. These 
analytical results can still be applied to unmatched regimes in the PWFA or in configurations 
associated with the external injection of electrons in plasma wakes. In fact, the beam energy, 
energy spread, emittance, transverse spot size, and number of accelerated electrons can only 
be effectively controlled when electrons are externally injected at the rear of the plasma 
wave [U [2S] • In this case it is advantageous to reduce the beam emittance to be as small 
as possible, thereby increasing the luminosity of the bunches, entering in a regime where 
Eq. ^ is valid. Later in this work we will show that this condition is also advantageous for 
beam polarization conditioning in plasma accelerators, and that the depolarization of lower 
emittance beams is lower in comparison to higher emittance bunches. 

In the conditions of Eq. ([T]), the electron trajectories are almost planar, and ds^/dt <^ 1 
which implies that s^ and s"^ + s^ = 1 — sL are conserved, with 5,^0 the initial tangential 



spin component. Under these conditions, the expression for the longitudinal component of 



the spin is found by combining Eqs. (4a) and (4b), yielding: 



Sz{t) = /l 



sin sin 



a -\ — ) Frdt + arctan I -^ 
V 7/ \s<f>o 



(8) 



where Szo, and s^o are the initial longitudinal and radial spin components, and where the 
radial component of the spin is given by s^(t) = 1 — s^q — s^(t). Eq. (p|, which is valid for 
arbitrary F^ as long as Eq. ([t]) is satisfied, describes the spin precession of a single electron 
in the PBA. 

The radial plasma focusing force in Eq. (Is]) depends on the electron trajectory as 
Fr = Fr{r{t)). The radial electron trajectory for low emittance electron beams ((|p±o|) ^ 
\/tt7(|x±o|)) is obtained from Eq. (|6|) by replacing x_L(t) by r(t), and x_lo by (xo|xo| + 



yo\yo\)/ y/xl + yl = tq. Thus Eq. (p| can be re- written as: 



Szilit)) = -^l-slosin 



ro$(7(t)) + arctan ' 

SrO 



(9) 



where $(7(t)) = -[(l+a7)(«^7o/7^)^/1 sin [2^/a{^ - v/7o)/^accci] • Equation Q gives the 
evolution of the longitudinal spin precession component for a single electron that propagates 
under linear focusing forces and constant accelerating gradients. It shows that Sz oscillates 
with the betatron frequency, and that the amplitude of the oscillations depends directly on 
the electron energy, since $ is a function of 7. Fig. [T] illustrates the evolution of Sz for a single 
electron according to Eq. ^ with ro = 0.5 c/up, 70 = 10^, Szo = 0.9, and s^o = 0.1, showing 
that the spin precession amplitude is on the order of s^o/100 thus suggesting that plasma 
accelerators can be efficiently used to accelerate polarized electron beams in conditions that 
are relevant for high energy physics experiments. 

The qualitative evolution of s{z) given by Eq. ^ is represented in Fig. 12} If the spin of 
the electron initially lies in the first or second quadrants of the [sr, Sz) phase space (i.e. if 



—IT < arctan (szo/sro) < 0), Sz reaches its maximum global extremum, Sz = -i/l — s^q, when 
7 is higher than 

[2 arctan (s^o/sro) - T^f ,-, r,\ 

7 > 7i = 1^ 4 4 2 ' \^^) 




FIG. 1: Evolution of the longitudinal spin component Sz as a function of 7 for a single electron 
with zero initial velocity, with tq = 0.5 c/ujp, 70 = 10^, Szq = 0.9, and s^o = 0.1. The solid lines 
denote the envelope of the oscillations. The dashed lines represent Sz{t). The amplitude of the 
oscillation is on the order of S20/IOO for an electron accelerating to 1 TeV, indicating the potential 
of plasma based accelerators to accelerate polarized electron beams. The period of the oscillation 
of Sz is equal to the betatron oscillations period. 



and it reaches its ininimuin global extremum, Sz = —\h- — s^n, when 7 is higher than 



7 > 72 



[2 arctan (s^o/sro) + ^r] 



(11) 



If the spin of the electron initially lies in the third or fourth quadrants of the (s^, Sz) phase- 
space, i.e. if TT > arctan (szo/^ro) > 27r, then 71 (72) is replaced by 72 (71). Inserting typical 



values for tq ~ c/ujp and for 70 ~ 10^ in Eqs. (10) and (11) yields 71 ~ 10^ and 72 ~ 10^"*^. For 
1 TeV electrons where 7 ~ 10^, this shows that the spin vector never performs one complete 
revolution in the (s^, Sz) phase space during the acceleration in PBAs. This feature is in 
contrast with circular accelerators, where the spin revolution frequency is of the same order 
of the orbital frequency of the electron beam [2] . 

The low amplitudes of oscillation for Szif) and Sr(t), are due to the low values that $(t) 
takes for typical PBA parameters. For $(t) ^ 1, Eq. (|9]) becomes: 



SM = S.0 + Sroro^t) - '^''om' ^ ^^^^^^3] 



(12) 



which reveals that s^ oscillates between the maximum s™^^ and the minimum s™'° given by: 



h {t) = ^zo ± s^o^oV " (77o) 7^ 



1/2 



(13) 



with a corresponding oscillation amplitude of: 



As^ ~ 2roS^oA/a(77o) a. 



\i/4, 



(14) 



Note that Eq. ( |14[ ) is valid as long as Szo or Sro remain in the same quadrant, i.e., 7 < 71. 



For 7 > 7i, it is modified to As!, = Asz + S, where 5 = ±a /l — s?q =F s™^'''™" (Fig. 



(a) 



(b) 



(c) 
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FIG. 2: Qualitative temporal evolution of Sz (solid-gray lines) during acceleration, (a) shows Sz{t) 



for for 7 < 71, (b) for 71 < 7 < 72, and (c) for 7 > 72- The solid-black lines are s 



max/min 



The 



circumferences show the possible trajectories of the spin in the {sz,Sr) phase space for constant 
S0. The arches on top represent the trajectory (s^, Sr) for the time frame correspondent to the left 
plots. 
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The evolution of As^ generally depends on the initial amplitudes of the betatron oscilla- 
tions (through To), on the reduction of the betatron oscillations amplitude as the electron 
accelerates (through 7), on the initial electron energy (through 70), and on the focusing 
force (through a). The spin precession amplitude Asz is lower for lower tq, i.e., for beam 
electrons that are closer to the axis. In addition, Asz increases throughout the accelera- 
tion, since the integral in Eq. (Is]) increases as the betatron amplitude of oscillation of single 
beam electrons is reduced during the acceleration j33]. If the betatron oscillations amplitude 
would be constant throughout the acceleration, the average value of As^ over a multiple of 
the betatron period would be constant since (As^) oc J^ia + 'y~^)Frdt = 0. Therefore, at 
the later stages of the acceleration, when the acceleration gradients are more significantly 
reduced, and when the variations of the amplitude of the betatron oscillations are thus more 
noticeable, d{Asz)/dt also decreases. Lower a's that correspond to weaker focusing forces 
also reduce As^. Therefore, the linear regime of plasma based acceleration can improve the 
polarization conditioning in plasma accelerators in comparison to the blowout regime. This 
requires, however, very narrow beams, such that the electron beam responds only to linear 
focusing forces. 

III. BEAM POLARIZATION 

The beam depolarization rates are critical to assess the potential of PBA for HEP ex- 
periments. The electron beam polarization component in a given spatial direction is given 
by the statistical average of s weighted by the polarized electron beam distribution function 

fern- 

P = f /esdxdvds / f /edxdvds, (15) 

We assume that the distribution function of the accelerated electron beam can be separated 
by its spatial (R(r)), velocity (V(v)), and spin (S(s)) distributions as: 

/, = R(r)V(v)S(s). (16) 

A transversely cylindrically symmetric Gaussian density distribution is considered, R(r) = 
nboexp[—2r^/a'^]Z{z), with n^ being the peak electron beam density, o",. the electron beam 



11 



transverse width, and where Z(z) is chosen to provide ideal beam loading [TO] which guar- 
antees that Ez can be the same throughout the beam. In fact, by choosing Z{z) according 
to [in], all beam electrons accelerate with the same accelerating gradient, and thus P is in- 
dependent of z. For a zero-emittance, and zero-energy spread electron beam, V(v) is given 
by V(v) = 5{'v±)5{vz — Vzo), where Vzo = (7o~l)^^^/7o is the initial beam longitudinal veloc- 
ity. The initial electron beam distribution function becomes fully characterized by choosing 
S(s); for the sake of simplicity, an uniformly polarized electron beam with S(s) = 6{s — So) 
is assumed. Beams with arbitrary polarization can be examined by allowing that So < 1, 
where |So| coincides with the initial polarization of the beam Pq. Physically this means that 
the average spin of the beam electrons located between r and r + 6r, and between v and 
V -|- (5v can be less than unity. Thus, the electron beam distribution function becomes: 



r 2r'^] 
fe = nbexp ^ Z{z)5(v±)5{vz-Vzo)5{s-So). (17) 

L ^r J 

The distribution function considered in Eq. ( |I7] ) which can be used as long as the beam 
emittance and energy spread is negligible, illustrates the key properties of spin precession in 
PBAs. 

The calculation of the beam polarization requires cartesian coordinates because the cylin- 
drical coordinate system is attached to the trajectory of each beam electron. Thus, in order 
to retrieve the beam polarization, {sr, s^j), Sz) must be projected in the cartesian coordinates 
basis, according to: 



s^ = Sr cos(0) - S0O sin(0) = ±Jl- s^- s^Q cos(0) - S0O sin (0) 



(18) 



Sy = s^o cos(0) + Sr sin {(f)) = s^o cos(0) ± W 1 - s2 - s^o sin (0), (19) 

being Sz identical in both coordinate systems. 

Analytical results for the beam polarization can be retrieved in the limit $ <^ 1, combin- 
ing Eq. Q, Eq. (|l5|), and Eqs.( [l7Hl9| , yielding: 



Pn 



El 
Po 



'0± 



1 r 



96 ^ ^ 



'S'oz 



$^2 <|.V,4 ^^,^. 

+ —r:r^ + 0{^- 



4 



48 



(20) 
(21) 
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where P_l = {Px,Py) is the polarization in each transverse direction, P^ the polarization 
along the z direction, So± = P±(t = 0), and Sq^ = Pz(t = 0). Comparison between 



Eqs. (20) and Eq. (21) reveals that |P_l| ~ {l/2)\Pz\. This indicates that the conditioning 
of transversely polarized electron beams with zero emittance is more effective than the 
polarization conditioning of longitudinally polarized electron beams. The total polarization 
of the beam P = a/P5_ + Pz is given by: 

+ "^^^^ ' '^^"" """"^^ +C>($'), (22) 



Po 8 384 



and confirms the aforementioned conclusion that the conditioning of transversely polarized 



zero-emittance electron beams is more effective since P/Pq is minimum for s^o = 0. Eq. (20) 
also reveals that the direction of P_l is preserved during the propagation, since Px/Py = 
Sxo/syo remains constant. The polarization vector then oscillates in the plane defined by 
s±o and e^. Similarly to the single electron spin precession dynamics, P never completes 
a revolution in the (P±,Pz) plane, and, in fact, the depolarization amplitude is lower than 
the spin precession amplitude, since |P| ~ s^ <^ 1. This reveals that the beam polarization 
can be maximized in the same manner that the spin precession amplitudes are minimized 
(cf. Sec.|lll). 

Fig. [3] illustrates the qualitative evolution of P, showing that the polarization oscillates 
with twice the betatron frequency between Pq and minimum values with an amplitude given 
by: 



Po 8 

|AP,| 5,0^2 («7o7)'/' a' 



Po 4 

I API {1 + si,) aUaiolY^' a' 



(24) 



(25) 



Po 8 

where AP_|_ = (AP^,, APj^) is the depolarization amplitude in the transverse directions. 



and AP is the total beam depolarization amplitude. Eqs. (23), (24), and (25) predict the 



maximum depolarization for an accelerated electron beam with a final energy E = meC^7 



at the end of the plasma accelerator. It is important to note that Eqs. (23), (24), and 



(25) imply that PBAs can be designed to keep the beam polarization exactly equal to the 

13 



initial polarization because for higher energy beams, with longer betatron wavelengths, it is 
conceivable that the length of the plasma can be adjusted to a multiple of half the betatron 
wavelength guaranteeing that P ~ Pq. 
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FIG. 3: Temporal evolution of the total polarization of a beam which accelerates from 0.5 GeV 
(70 = 10^) to 1 TeV (70 = 2 X 10^), considering Pzo/Pq = 0.9. The solid black lines are the upper 
and lower envelopes AP/Pq. The dashed lines represent the evolution of AP/Pq. 



IV. SELF-GUIDED VS. EXTERNALLY-GUIDED LWFAS 



Two regimes for stable propagation have been proposed in the blowout regime of LWFA, 
either using external guiding structures (such as parabolic plasma channels), or by adjusting 
the laser and plasma parameters to ensure self-guiding [9]. It is straightforward to retrieve 
the total depolarization associated with externally-guided or self-guided LWFAs of similar 
output energy according to the scaling laws from Ref . [9] . It was assumed in this calculation 
that the physical dimensions of the beam are the same for the externally-guided or the self- 
guided scenarios, while the rest of the parameters such as laser energy, and plasma density, 
were chosen to ensure stable accelerating structures. Denoting the total depolarization 
associated with the self-guided propagation regime by AP^^/Pq, and the depolarization 



associated with the externally-guided propagation regime by AP'^^/Pq, Eq. (25) leads to: 
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— — ~ 1.17 X 10"^^ f' ^ '" '— -^ , (26) 

for the self-guided propagation regime, and as 

^^''.,8.9xlO-<i±%^f,^Vf^V. (27) 



Pa ' 7"= VlOfim/ V Ao / 

for the externally-guided propagation regime, where Aq is the laser central wavelength in 
/xm, and a^ is the electron beam width measured in /im. Although the depolarization 



increases during the acceleration, Eqs. (26) and (27) reveal that higher energy matched 
propagation LWFAs lead to lower depolarizations. This is because a higher energy matched 
LWFA uses lower densities, lowering the normalized 0"^ (the beam dimensions in physical 
units are kept constant), and thus AP/Pq. The focusing force is also lower for externally 
guided propagation regimes, which further decreases the depolarization rates. In standard 
accelerators, however, AP/Pq increases for higher energies. This result thus emphasizes an 
important advantage to accelerate highly polarized electron beams to very high energies in 
PBA. 



The results of Eqs. (26) and (27) can be illustrated by considering typical electron beam 
parameters. At SLAG, &,. = 10 fim, and 70 = 6 x 10^. Thus, the maximum total depolar- 
ization associated with a 0.5 TeV accelerator is AP'^^/Pq = 0.16% and AP^s/Pq = 0.025%, 
where we used a = 0.35, estimated from 3D PIC simulations [M]. We note that these re- 
sults do not take into account the beam depolarization due to the beam-beam interaction at 
the interaction point. Nevertheless, it is instructive to compare the typical depolarizations 
expected in PBAs with the expected depolarization at the ILC, which is on the order of 
^piLC ^ Q ]^ _ 2 1^ [-we note that AP^^'~' is mostly due to the depolarization at the beam- 
beam collision point]. Although this work only takes into account the beam depolarization 
during the acceleration, this comparison thus reveals that plasma accelerators should also 
be able to fulfill the depolarization requirements of HEP experiments. 

V. NUMERICAL SIMULATIONS 

In order to investigate the acceleration of polarized electron beams in realistic conditions, 
the T-BMT equations are further explored with a numerical spin-tracking algorithm. The 
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evolution of the polarization is firstly examined using prescribed fields and electron trajec- 
tories (cf. Sec. In]), where very fast parameter scans can be performed even for > 1 TeV 
accelerators. Then, the results from one-to-one 3D QuickPIC [32] simulations are presented, 
providing strong evidences for PBA with low depolarization, in agreement with the theoret- 
ical predictions. 

The similarities between Eq. ^ and the magnetic field component of the Lorentz force 
suggest that a non-relativistic Boris-pusher |35] can be used to solve the T-BMT equations. 
The rotation of the spin is determined by f2, in the same way the rotation of the electron 
velocity is determined by B in the Lorentz force equation. For At/A/3 ^ 1, the spin rotates 
by an angle 6, given implicitly by d = 17 tan (5/2) ~ f2At/2. The rotation of s between 
each time step is given by: 



Sj + Sj X d. 



(28) 



and 



'i+l 



Si + 



2s\ X d 

l + |d|2 



where i is the iteration number. Eqs. (28) and (29) guarantee that |s 



Sj+l 



(29) 
P, providing 



an accurate numerical resolution of the T-BMT equations, to second order in At. 



In order to compare the analytical model with Eqs. (28) and (29), we consider prescribed 



radial focusing forces, F^ = ax^, and the electron trajectories described by Eq. ([6]). It then 
follows that 



and 



n. 



— a 



lit) 



—a a + 



lit) 



(Ex — By) 

x{t), 



ny = (a + -iy^ (Ey - B,) 

y{t). 



a \ a + 



lit) 



Fig. |4| shows a comparison between 5™°/™'^'' (^gee Eq. (13)), AP/Pq (see Eq. (22)), and 



the numerical resolution of the T-BMT Equations (Eqs. (28) and (29)). A zero-emittance 
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electron beam was initialized with ar = 0.1 c/up and 70 = 10^, which accelerates with 
-^accei = 0.6 during L = 2.4 x lO'* c/up. The numerical results are in excellent agreement 
with the analytical predictions. 
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FIG. 4: Comparison between the analytical model and the numerical resolution of the T-BMT 
equations using prescribed fields, (a) Longitudinal component of the spin evolution for a beam 
with ar = 0.1, which accelerates in a constant accelerating gradient -Eaccei = 0-6, with a = 0.35. 
(b) Evolution of the total depolarization of the beam. The solid line corresponds to the analytical 
result, while the dashed line corresponds to the numerical solution. 



Using the spin-tracking numerical algorithm, the depolarization associated with realistic 
electron beams, with finite emittance, can be investigated. Of particular relevance is the 
assessment of the role of the beam emittance for typical electron beams now available. 
Fig. [5] shows the evolution of the polarization of a SLAC-like electron beam, initialized with 
ar = 10 /im and 70 = 6 x 10*^. The beam normalized emittance is eNx = 50 mm ■ mrad and 
ejsfy = 20 mm-mrad in the x and y directions, reaching 0.5 TeV at the end of the acceleration. 
The plasma parameters were chosen to mimic LWFA stable propagation regimes: for the 

15 021-3^ -Eaccei = 0.7, and a = 0.35. The longitudinal 
depolarization rates are only marginally increased due to the non-zero emittance, which 
plays a stronger role in the external-guiding scenario. However, the total depolarization rates 
increase due to a drift of the transverse spin component of single beam electrons, which leads 



self-guided case uq = 1.25 x 10^^ cm ^, -Eaccei = 1.65, and a = 0.5; for the externally- 
guided scenario, uq = 2.4 x 10^^ cm~^, -Eaccei 



to stronger variations of P_l (cf. Eq. (23)). The numerical results also show that the effects 
of the beam emittance on the depolarization can be neglected for the early propagation. 
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The expected total depolarization associated with the ILC, including all the depolarization 
sources present at a linear collider, is also shown in Fig. [5] The main depolarization source 
at the ILC conies from the beam-beam interaction at the interaction point, which provides 
typical depolarizations on the order of 0.1 — 0.2%. Although this work does not include the 
depolarization that occurs at the collision of the two colliding beams. Fig. [5] suggests that 
the main source of depolarization in a plasma based collider may also come from the beam- 
beam interaction at the interaction point provided that the beam emittance is sufficiently 
low, so that the beam depolarization is kept below 0.1 — 0.2% during the acceleration. 

3D PIC simulations using the quasi-static PIC code QuickPIC [32] have been performed 
to examine in detail the spin precession and polarization dynamics in LWFA scenarios. 
To this purpose, a random sample of the externally injected electrons was tracked using a 
particle-tracking algorithm [36], which stored the trajectories and fields associated with the 
accelerated electron beam. 

The simulations, which worked in the externally guided propagation regime, used the 
matched parameters from Refs. [HI HO] for a 800 J laser pulse. The laser central wavelength 
was Ao = 0.8 /im, normalized peak vector potential oq = (eAo/meC^) = 2, spot size Wq = 
140 /im, and duration tpwhm = 311 fs. The transverse plasma profile is given by n{r) = 
no[l + An{r/WOy] for r < re, falling linearly for r > Tc, where hq = 1.15 x 10^^ cm~^. 
An = 0.375, Tc = 560 yum, and Vf = 2380 fun. The size of the simulation box is 2380 /im x 
2380 /im X 450 /xm, divided into 128 x 128 x 1024 cells, with 4 particles per cell. The 
externally injected electron beam was initialized at the back of the first plasma wave, in the 
region of maximum accelerating fields, with zero-emmitance, and zero-energy spread. The 
beam charge was g = 80 pC, and 70 = 10'^. The beam density is transversely Gaussian with 
width o"r- = 23 /im and length a^ = 23 /im. In the simulations, 125 electrons were followed 



in time, and the particle-tracked data was post-processed according to Eqs. (28) and (29), 
in order to retrieve the evolution of the beam polarization. 

Fig. [6] shows a comparison between the simulation and the analytical model. For the 
analytical model results, -Eaccei = ^^0/2 ~ 0.7 [9] was used, and the coefficient a = 0.35 
was retrieved from the simulations (note that even though in the simulation a parabolic 
plasma channel is used, the focusing force is still linear in the region of the beam). Excellent 
agreement is obtained for s^ until 3 m of propagation (Fig. M^^a)), when 7 ^ 5 x 10'^ (~ 
25 GeV). Despite the not perfectly uniform acceleration along the electron beam as the 
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FIG. 5: Beam polarization evolution for a 0.5 TeV PBA. (a) Evolution of the longitudinal com- 
ponent of the polarization in a self-guided regime, (b) Evolution of the total polarization in a 
self-guided propagation regime, (c) Evolution of the longitudinal component of the polarization 
in an externally guided propagation regime, and (d) Evolution of the total polarization in an ex- 
ternally guided propagation regime. The gray (black) line corresponds to a SLAC-like beam with 
(without) emittance, and the black dashed line corresponds to the typical depolarization that is ex- 
pected for the ILC, which includes the depolarization of the beam at the interaction point between 
the colliding beams. 



acceleration progresses (-Eaccei drops), very good agreement was also found for the evolution 
of the polarization (Fig. M^b)). We note that the oscillation amplitude of the polarization 
does not match the analytical model because of the limited number of tracked particles in 
the 3D simulations, which is limited by computational constraints. However, the trend of 
the polarization is fully captured. For propagation distances larger than 3 m the model 
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overestimates the spin precession amplitude because of the reduction of the acceleration 
gradient (due to the laser pump depletion and dephasing), and also because the focusing 
force changes (a decreases during the simulation). 
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FIG. 6: QuickPIC simulation results for a externally guided LWFA, using a 800 J laser pulse to 
accelerate a zero-emittance electron beam to 70 GeV. The predictions of the model are shown in 
black, and the results of the simulation shown in gray, (a) shows the evolution of Sz for a randomly 
chosen beam electron, (b) shows the evolution of the total polarization. The black curves are the 
analytical results for the maximum and minimum values of AP/Pq. (c) shows the evolution of the 
electron beam energy, (d) shows the radial trajectory of the same electron of (a). 



In order to identify the role of the beam emittance in more realistic scenarios, Quick- 
PIC [32] simulations using non-zero electron beam emittances were performed, using e™ = 
62.5 mm-mrad (typical SLAC beam emittance) and keeping the rest of the laser, plasma, and 
beam parameters unchanged. For these values, Eq. (H) yields e^^^ = 167 mm-mrad ^ 3e^™. 
Fig. [TFa) shows the evolution of s^ for a given beam electron. The initial electron velocity 
shifts the value of Sz from the theoretical prediction because the amplitude of the electron 
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oscillation changes due to the initial electron velocity. The total beam polarization, shown 
in Fig. [Wb), indicates that the depolarization rates are slightly higher in comparison to those 
of Fig. [61 
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FIG. 7: QuickPIC simulation results for a externally guided LWFA, using a 800 J laser pulse to 
accelerate an electron beam with normalized emittance e^ = 62.5 mm- mrad to 70 GeV. The black 
curves refer to the predictions of the analytical model, and the gray curves represent the simulation 
results (a) shows the evolution of Sz for a given beam electron. The results of the model assumed 
no initial transverse electron momentum, (b) shows the evolution of the total polarization. The 
analytical results illustrate the maximum and minimum values of AP/Pq assuming a zero-emittance 
electron beam. 



VI. CONCLUSIONS 

This Paper describes the key features associated with the acceleration of polarized electron 
beams in PBAs. It was found that the spin precession dynamics in the PBA is closely 
related to the betatron oscillations in the ion channel, and that for zero emittance beams, 
the polarization vector oscillates in the plane defined by the propagation direction and the 
initial polarization vector direction. Therefore, if the beam is initially polarized in an optimal 
direction for a certain HEP application, it remains in that direction during the acceleration 
stage. Furthermore, higher energy matched LWFAs lead to lower depolarizations. 

Our model was validated by the numerical resolution of the T-BMT equations, using 
prescribed electromagnetic plasma wave fields and electron trajectories, and by using the 
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fields and electron trajectories retrieved from 3D PIC simulations in QuickPIC [32]. These 
results further confirmed the potential of PBA to accelerate polarized electron beams to high 
energies with final depolarizations which can fulfill the requirements of high energy physics 
experiments provided that low beam emittances are used. 

The analytical work presented in this Paper was obtained in the limit of small emittances 
such that ejv[mm ■ mrad] < e^'"^'^ = 188v^(a^/10 /im)^ (^o/lO GeV)^/^ K/IO^^ cm-^)^/^ 
We note that the beam luminosities are also larger in the conditions where ejv[mm-mrad] ^ 
gmrad Therefore our analytical model describes the spin precession in PBAs in the conditions 
that are directly relevant for HEP applications. Moreover, numerical modeling showed that 
larger initial beam emittances lead to larger depolarization rates, specially in the transverse 
directions. 3D PIC simulations in QuickPIC also indicated that the depolarizing effects due 
to the beam emittance are balanced by the decrease of the accelerating gradients caused by 
the laser pump depletion. 

In future works it could be important to further assess the effects of the beam emittance, 
beam energy spread in the beam depolarization, along with the depolarization associated 
with the beam-beam interaction at the interaction point in plasma based accelerator collid- 
ers. 
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